Persistence of coherent quantum dynamics at strong dissipation 



o 

(N 



o 

CD 



O 
> 

ON 

O 

o 



X 



Denis Kast and Joachim Ankerhold 
Institut fur Theoretische Physik, Universitdt Ulm, Albert- Emstein-Allee 11, 89069 Ulm, Germany 

(Dated: January 10, 2013) 

The quantum dynamics of a two state system coupled to a bosonic reservoir with sub-Ohmic 
spectral density is investigated for strong friction. Numerically exact path integral Monte Carlo 
methods reveal that a changeover from coherent to incoherent relaxation does not occur for a broad 
class of spectral distributions. In non-equilibrium coherences associated with substantial system- 
reservoir entanglement exist even when strong dissipation forces the thermodynamic state of the 
system to behave almost classical. This may be of relevance for current experiments with nanoscale 
devices. 
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Introduction-The impact of dissipative environments 
on the dynamics of quantum systems has regained inter- 
est due to the boost of activities to tailor atomic, molec- 
ular and solid state structures with growing complex- 
ity [Il-Q- Common wisdom is that quantum coherence 
is inevitably destroyed at sufficiently strong coupling to 
broadband reservoirs. A paradigmatic model is a two- 
state system interacting with a reservoir of bosonic de- 
grees of freedom (spin-boson model) whichplays a funda- 
mental role in a variety of applications [J-l6| . At low tem- 
peratures and weak coupling, an initial non-equilibrium 
state evolves via damped coherent oscillations, while at 
stronger dissipation it displays an incoherent decay to- 
wards thermal equilibrium. This is often understood as 
a quantum to classical changeover. The question we pose 
here is whether this picture always applies. 

The spin-boson (SB) model has recently gained re- 
newed attention for reservoirs with sub-Ohmic mode dis- 
tributions Js (uj) oc a where a denotes a coupling con- 
stant and < s < 1. This class of reservoirs constitutes 
a dominant noise source in solid state device^ at low tem- 
peratures such as superconducting qubits [T] and quan- 
tum dots 8] with the spectral exponent s determined by 
the microscopic nature of environmental degrees of free- 
dom. It also appears in the context of ultra slow glass dy- 
namics , quantum impurity systems jlQ] , and nanome- 
chanical oscillators [ll|. Advanced numerical techniques 
12 - Tsi] have revealed that at zero temperature the equi- 
librium state exhibits at a critical coupling strength ac 
a quantum phase transition (QPT) from a delocalized 
phase with tunneling between the two spin orientations 
(weak friction) to a localized one with almost classical 
behavior (stronger friction) [l^ . 

The time evolution of the polarization towards 
these asymptotic phases shows coherent oscillations or 
classical-like monotone decay. It has been argued that 
with increasing friction the relaxation dynamics always 
turns from coherent to incoherent fl3|. Indeed, numeri- 
cal studies [Til confirmed this picture for reservoirs with 
1/2 < s < 1, but the situation for < s < 1/2 remained 
unclear, mainly because approaches used previously are 



restricted to the regime of weak to moderate dissipation 
12, [l^. The goal of this Letter is to attack this lat- 
ter regime via real-time path integral Monte Carlo tech- 



niques (PIMC) [20j, |2l| which cover also strong friction. 



We verify that coherences persist in non-equilibrium for 
arbitrary coupling strength to a heat bath even when the 
thermal state is essentially classical. Our results shed 
new light on our understanding of the quantum-classical 
crossover and may be accessible experimentally. 

Reduced dynamics - We consider a symmetric SB 
model 



hA 



(1) 



with a two state system (TSS) which interacts bi- 
linearly with a harmonic reservoir via the bath 
force £ = J2a^a{ba + AH relevant observ- 

ables are obtained from the reduced density operator 
p{t) = TrB{exp{-iHsBt/h)W{0)exp{iHsBt/h)}, where 
the initial state has the form 

M/(0)=ps(0)e-^(«--^)/ZB. (2) 

Here, according to typical experimental situations Q the 
bath distribution is equilibrated to the initial state of 
the TSS which is ps{0) = \+l) (the eigenstates of 
az are |±1) with the bare tunneling amplitude A between 
them); the partition function is denoted by Zb- We are 
interested in the real-time dynamics of the observables 

P^(t) = {(Ju{t)) = Tr {a^ p{t)} , iy = x,y,z, 

where Pz describes the population difference (polariza- 
tion) and Px the coherence between the sites | ± 1). 

A non-perturbative treatment is obtained within the 
path integral formulation. The Pi, is expressed along a 
Keldysh contour with forward a and backward tr' paths 
5] . The impact of the environment appears as an influ- 
ence functional introducing arbitrarily long-ranged inter- 
actions between the paths. Switching to 77/^ = cr ± tr' 
one arrives at 
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FIG. 1: Domains of coherent (white) and incoherent (shaded) 
dynamics of the SB model for a sub-Ohmic environment Q 
at zero temperature. Above (below) the solid line ac{s) the 
system asymptotically reaches a thermal equilibrium which is 
localized (delocalized) [23|. A coherent-incoherent changeover 
only occurs along the dashed line aci{s) (0. Dotted lines 
refer to results in Figs. (2] [S] respectively. 



with the contribution Ai, in absence of dissipation and 
the influence functional 

^ivA]^ [duf dvi{u)[Q'{u~v)i{v) + iQ"{u-v)ri{v)]. 
Jo Jo 

According to ^ one sums over all paths with |r](0)| — 
1,^(0) = and T]{t) = l,C{t) = for and r]{t) = 
0, |^(<)| ^ 1 for P^, respectively. The kernel Q = Q' + iQ" 
is related to the bath correlation Q{t) = {£{t)£{0)) /h^ 
and is thus completely determined by the spectral func- 
tion J{uj) = nJ2a ^a^^ ^ ^a)- Foi sub-Ohmic spectral 
distributions 



< s < 1, 



(4) 



Jsiuj) = 27raa;,i-'^w''e-'^/'^= 
one has at zero temperature Q 

Qo{t) = 2ar(s - 1)[1 - (1 + iuj.ty-'] 



In the Ohmic case s ^ I the coupling a coincides with 
the Kondo parameter K. The cut-off ujc corresponds e.g. 
to a Debye frequency [5f or to a parameter of an electro- 
magnetic environment [8|. 

A direct evaluation of is extremely challenging due 
to the retardation in the influence functional which grows 
with decreasing temperature. In this situation PIMC 
methods are very powerful means to explore the non- 
perturbative range including strong coupling [2l| . 
Here, we extend this formulation to simulate not only 
populations Pz but also coherences Px (see (23|). 

Population dynamics- The dynamics of Pz (t) directly 
displays the impact of decoherence as the bare two 
state system is completely quantum and has no classi- 
cal limit. The expectation is that finite friction induces 
damped oscillations on a transient time scale of the form 
Pz{t) ^ e~''^ cos{nt) . As long as 7^ the system dy- 
namics is said to be coherent, otherwise it is incoherent. 




FIG. 2: Dynamics of Pz{t) at T = and A/loc = 0.1 according 
to PIMC for s — 0.75 and couplings a with Qc ~ 0.3 > aci ~ 
0.22 (purple line). Blue (red) lines refer to coherent a < aci 
(incoherent a > aci) dynamics (cf. Fig. [1]). Statistical errors 
are of the size of the symbols. 



For vanishing friction Pz{t) — cos(At), while finite cou- 
pling leads asymptotically to a thermal equilibrium with 
Pz{t 00) ^ (delocalized) or Pz{t ^ 00) ^ Q (local- 
ized) as shown previously [l3| . 

To analyze the dynamical features in detail (cf. Fig.[T|), 
we consider the case of zero temperature and cjc 3> A, 
and start in Fig. [5] with Pz{t) in the range 1/2 < s < 1 of 
spectral exponents. For s fixed, PIMC simulations dis- 
play with increasing coupling a a changeover from coher- 
ent to incoherent motion, i.e., from damped oscillations 
to overdamped decay |18| . It appears in that domain 
in parameter space where the thermal state is delocal- 
ized (Fig. [T]) . Accordingly, the critical coupling strength 
for the coherent- incoherent turnover aci{s) is always 
smaller than the critical coupling ckc(s) for the QPT. In 
limiting cases one finds ads — 0.5) w aciis = 0.5), 
while for an Ohmic bath (s = 1) the known result is 
confirmed a^s = !)«!> aci{s = 1) w 0.5. 

In contrast, simulations in the range < s < 1/2 reveal 
a different scenario. Oscillatory patterns survive even for 
coupling strengths far beyond the critical coupling ac{s) 
for the QPT (see Figs. [II E]). PIMC data up to ultra- 
strong couplings a — 30ac do not show a changeover to a 
classical-like decay for exponents up to s = 0.49. The os- 
cillation frequencies fls{cx) in Pz{t) increase with increas- 
ing coupling and exhibit a scale invariance according to 
ns(/ca) — K^lsia) (inset Fig. [3] and [l^l). Two observa- 
tions are intriguing: (i) even in the regime where friction 
is so stron g th at the thermal state resides in the classical 
phase IJ, |l5| (spin almost frozen with Pz 



1 > Px 



does the non-equilibrium dynamics of the TSS preserve 
quantum coherence and (ii) the domain where these co- 
herences survive, covers a broad range of spectral distri- 
butions up to exponents s close to s = 0.5 (see Fig.[T]). 

To gain analytical insight, an approximate treatment 
is provided by the non-interacting blip approximation 
(NIB A) [5]. The starting point are exact equations of 
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FIG. 3: Same as in Fig.[2]but for s — 0.25 and values a > 
0.022. The inset shows a blow-up for very large couplings. 

motion that can be derived from ([3]) p|, i.e., 



and decay rate 71 (a) = aA tt ■\/wc/4A. There is no 
pole with A = meaning that Pz{t 00) — >■ (de- 
localization). Apparently, the oscillation frequency de- 
creases with growing coupling a while the damping rate 
increases. This indicates a changeover from coherent mo- 
tion with J7i /71 > 1 to incoherent motion with ilj. = 
2 ' 2 2 

for larger a. 

(ii) In the domain s ^ 1 one may use the expan- 
sion Qo{t) w Astexp[iTT{l — s)/2] with coefficient = 
dLuJsiuj)/{TTLj) = 2awc|r(s- 1)|. Then, P(A) has one 
pole at A = so that Pz{t) relaxes towards an asymp- 
totic limit Pz{t 00) > (localization). The other two 
poles are complex conjugate where the imaginary part 
describes oscillations in (t) with frequency 



ns{a) w AsSin[(l - s)tt/2] 



2aujr 



(8) 



Pzit) = - / duK,{t-u)Pziu). 
Jo 



(5) 



Within NIBA the kernel is evaluated as ~ Km , where 



ifjv,.(t) = AV3'(*) cos[Q"(<)] 



(6) 



Let us now analyze how incoherent decay for Pz{t) may 
appear out of ([5]). In KN,z{t) the correlation Qgit) 
induces damping while Qo(i) is responsible for oscilla- 
tory motion. Both functions are positive and monoton- 
ically increase in time leading to damped oscillations in 
K^^zit)- If this damping is sufficiently strong on the time 
scale of the bare dynamics 1/A, i.e. if the kernel is suffi- 
ciently short-ranged in time, the population Pz{t) decays 
incoherently. We estimate this to be the case if at the 
first zero t = t^, oi the kernel where Qo(i») = 7r/2, the 
damping obeys (5o(i*) > 1. One finds this condition to 
be fulfilled only for s > 1/2. There, the borderline be- 
tween the two dynamical regimes is determined by that 
coupling strength at which = 1/A, i.e.. 



aciis) 



A 



4|r(s - l)|cos(s7r/2) 



(7) 



This expression includes the known Ohmic result 
0^0/(5 = 1) = 0.5 and captures accurately numerical data 
from our PIMC. Further, it confirms our numerical find- 
ing that ac{s) > aci{s) (see Fig. [T]). For fixed s > 0.5, 
increasing friction first destroys coherent dynamics before 
it asymptotically induces also localization. In contrast, 
for spectral functions < s < 0.5 a coherent- incoherent 
changeover never appears. 

In principle, the oscillation frequency ils and the 
damping 7^ of Pz (t) can be extracted as complex- valued 
poles from the Laplace-transform of ([5]), i.e., P(A) — 
1/[A + Kn,zW]- However, specific results can be found 
only for (i) s — 0.5 and (ii) s ^ 1 (for details see [23(). 

(i) In the regime a <^ y/ A/wc, the poles are given 
by A^ = iifii — 71 with frequency Qi{a) « A — 71 



The real part corresponds to a damping rate 7^ (a) « s A^ 
which saturates for s — > at 70 = 2aujc- Hence, for 
fixed s, both rate and frequency increase with increasing 
a such that the dynamics of the TSS is always under- 
damped fls/js — 1/s^ 1- Quantum coherent dynamics 
persists up to arbitrarily large couplings in accordance 
with the numerical data (cf. Fig. [3|). Further, the above 
result confirms the scaling property of fig (a) found al- 
ready in the PIMC. Note that details of the cut-off pro- 
cedure may enter via A^ only. Since it is determined 
by the low frequency behavior of J(u;), the dynamics is 
independent of the cut-off scheme for Wc 3> A as also 
seen in the PIMC (not shown). Similarly, the bare fre- 
quency of the TSS has disappeared. It only governs the 
dynamics for ultra-short times t <^ 1/wc before reser- 
voir modes can respond. There, the bare Schrodinger 
dynamics predicts the universal quadratic time depen- 
dence Pz{t) « 1 — A^t^/2. After this initial delocaliza- 
tion process, the reservoir tends to prevail with the low 
frequency modes in £ [see ([T])] acting on the TSS effec- 
tively as a static energy bias of strength —hAs according 
to the coupling —azS-12. This is in agreement with the 
observation that TiAg corresponds to the energy needed 
to re-organize the reservoir once the spin has flipped 24 1 . 

Coherences and entanglement- To further verify the 
quantum nature of the oscillatory population dynamics, 
we also monitor Px{t) (Fig. |4]). Indeed, one observes 
substantial quantum nonlocality which also suggests sub- 
stantial entanglement between the TSS and its surround- 
ings. At T = this entanglement can be extracted 
from the von Neumann entropy S{t)/kB — — w+ln(w+) — 
ii;_ln(ii;_) with w± = i[l ± y'P4t)^ + Py{t)^ + Pz{tf] 
(cf. Fig. 0]). PIMC simulations reveal strong entangle- 
ment well in the regime a > Oc (localization) on time 
scales where oscillatory patterns in Pz{t) and Px{t) occur. 
For longer times and stronger coupling, entanglement de- 
cays montonically. This is in agreement with previous 
findings in thermal equilibrium: entanglement tends to 
zero in the localized phase for couplings somewhat above 
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FIG. 4: Coherences (left axis) and entropy S{t) (right 

axis) for s = 0.25 (a^ ~ 0.022) at T = according to PIMC. 




Initial preparation- Persistence of coherence is associ- 
ated with a strong dependence on initial preparations. 
Of particular interest for experimental realizations are 
preparations, where the reservoir is out of equilibrium 
with respect to the initial state of the TSS. This means 
that in ^ one must replace 8 fiE with /.j ^ 1. The 
cloud of bath modes is shifted to the left (right) for /i < 1 
(/i > 1) with /i = being the equilibrium orientation of 
the bare bath. According to Fig. [5l with decreasing /i 
the effective oscillations frequency of the TSS decreases 
while at the same time the initial loss in population in- 
creases. This behavior can be understood from the fact 
that this preparation can equivalently be described by a 
time dependent bias e^(t) = (1 — fi)Q"(t) of the TSS, i.e. 
by an additional term he^{t)az/2 in ([ij 25]. In particu- 
lar, for reservoirs with spectral exponents s ^ 1 one has 
Q"{t) « As so that ~ (1 — fj,)As becomes static. Effec- 
tively, this bias adds to the bias induced by the sluggish 
modes in the system-reservoir coupling —£az/2 [see be- 
low ([8|)] to produce a net bias etot ~ — M^s- As a result 
the TSS regains bare dynamical features for /i — )• 0. This 
analysis can now be extended to other preparations of 
the TSS alone (superpositions of | ± 1)) which, however, 
display qualitatively a similar picture. The same is true 
for asymmetric TSS and/or finite temperatures as long 
as corresponding energy scales do not exceed the bare 
tunneling amplitude. Results will be shown elsewhere. 

Summary-We have shown that in non-equilibrium co- 
herent dynamics can persist for ultra-strong coupling to 
a broadband reservoir. For a SB model in the sub-Ohmic 
regime, stronger friction does not induce incoherent re- 
laxation for spectral exponents < s < 1/2 even when 
the thermal equilibrium is almost classical. These find- 
ings shed light on our understanding of decoherence in 
open quantum systems and are thus of relevance for 
current experiments in nanoscale structures. The case 
s = 1/2 can be realized through a charge qubit sub- 
ject to electromagnetic noise [8|, while recent progress 



FIG. 5: Pz{t) at s — 0.25, q = 0.1 for various initial prepara- 
tions of the reservoir with ujc/A — 10. See text for details. 

in engineering local environments for Cooper pair boxes 
may allow to study cases with s < 1/2 [26[. Reservoirs 
with s — ^ have been proposed as models for 1/f noise 
in superconducting circuits where temperature enters as 
an effective parameter Q. While in this latter case de- 
tails need to be investigated, a promising alternative are 
trapped ion systems as shown in p7| . 
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